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Abstract 

It has been previously observed that during a pumping test in heterogeneous me- 
dia, drawdown data from different time periods collected at a single location produce 
different estimates of aquifer properties and that Theis type-curve inferences are more 
variable than late-time Cooper- Jacob inferences. In order to obtain estimates of aquifer 
properties from highly transient drawdown data using the Theis solution, it is neces- 
sary to account for this behavior. We present an approach that utilizes an exponen- 
tial functional form to represent Theis parameter behavior resulting from the spatial 
propagation of a cone of depression. This approach allows the use of transient data 
consisting of early-time drawdown data to obtain late-time convergent Theis parame- 
ters consistent with Cooper-Jacob method inferences. We demonstrate the approach 
on a multi-year dataset consisting of multi-well transient water-level observations due 
to transient multi-well water-supply pumping. Based on previous research, transmis- 
sivities associated with each of the pumping wells are required to converge to a single 
value, while storativities are allowed to converge to distinct values. 



1 Introduction 



Aquifer property inferences obtained using the Theis type-curve method ( Jacob , 1940) (Theis 



method) and the Cooper- Jacob straight-line approximation method ( Cooper and Jacob 



1946) (Cooper- Jacob method) at a given location have been observed to differ (Ramey 



1982 



Butler, 1990). Theoretical investigations by Dagan (1982) utilizing a perturbation 



expansion approach on idealized scenarios demonstrate that effective hydraulic conductivity 
(transmissivity in 2D) decreases from the arithmetic mean conductivity to a convergent value 
over time. More recent numerical and field investigations demonstrate that Theis solution 



parameters ( Theis , 1935 ) estimated at a location at various times during a pumping test have 



been observed to decrease at early times converging to stable values at late-times ( Wu et al. 



2005 Straface et al, 2007). Butler (1990) contributes this characteristic of Theis solution 
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parameters to the fact that at early times, while the cone of depression is approaching the 
observation location, the drawdown is affected by many factors, such as: skin effects; well 
loses; and aquifer heterogeneities encountered by the cone of depression, complicating the 
estimation of stable parameters. However, at late times when quasi-steady state conditions 
have developed (i.e. when pressure gradients have reached steady state but pressures remain 
transient), the stable parameter estimates are consistent with aquifer property inferences 
that would be obtained using the Cooper- Jacob method. This implies that the late-time 



parameter estimates provide interpreted aquifer properties (as defined by Sanchez- Vila et al. 



(2006)) representative of the support scale defined by the distance between the pumping and 



monitoring wells (Neuman, 1990 Neuman and Di Federico, 2003). 



Obtaining variable model parameter inferences indicate the inadequacy of a model to 
represent a system, as parameters are intended to represent invariant intrinsic properties 
of the system. The limitations of applying the Theis solution to model typical pumping 
tests is not a matter of debate, as its inadequacies are readily apparent by the assumptions 
required in its derivation (Theis, 1935) (e.g. fully penetrating well, infinite lateral extents, 



homogeneous properties, unperturbed initial conditions, confined aquifer). Recognizing these 
limitations, the question becomes whether or not the model can be useful. We agree with 
previous researchers that the Theis solution is useful for obtaining aquifer property inferences 
that characterize the groundwater transport if late-time drawdown data is used consistent 



with the Cooper- Jacob method ( 


Butler , 


1990 


Meier et al. 


1998 


Sanchez- Vila et al. 


1999 


Knudby and Carrera, 


2006 


Trinchero et al. 


c 


2008 


). As noted by 


Butler 


(1990 


) in reference 



to the use of the Cooper- Jacob method, the advantage of drawing inferences from late-time 
drawdown data is that the estimated parameters will be independent of the numerous early- 
time effects that can influence the drawdown at the initial stages of expansion of the cone of 
depression. 

Obtaining late-time pumping data at quasi-steady state is not always possible, however, 
as it may not be feasible to continue a pumping test for a sufficient duration to allow 
quasi-steady state conditions to develop. Or, in cases where an existing water-supply and 
water-level elevation dataset is available from a municipal water supply network, quasi-steady 
state may not be reached due to a high frequency of cycling multiple supply wells on and 
off to: meet shifting demand; to take advantage of lower cost off-peak electrical rates; and 
perform well maintenance and/or repair. In this paper, we present an approach that allows 
convergent parameters to be obtained from transient drawdown data by accounting for the 
behavior of Theis parameters at early times. 

The proposed approach is demonstrated on a long-term highly-transient drawdown record 
collected at the Los Alamos National Laboratory (LANL) site where the water-level tran- 
sients result from multi-well municipal water-supply pumping. The pumping regimes are 
highly transient, cycling diurnally and seasonally, including variations due to maintenance, 
repair, and shifting supply loads within the network. As a result, the drawdown at monitor- 
ing wells within the network do not fully attain quasi-steady state as new pressure influences 
begin to propagate through the aquifer as the pumping wells cycle on and off (Harp and 



Vesselinov 2010). The use of a long-term dataset containing multiple pressure influence 
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cycles has certain advantages, such as: reduction of measurement errors; improved charac- 
terization of the hydraulic response allowing the refinement of hydrogeologic inferences; and 
the lack of the expense and coordination of a conventional pumping test. We demonstrate 
the inference of aquifer properties from this dataset by considering the transient early-time 
behavior of Theis solution parameters. 

As the approach presented here utilizes observations, numerical experiments, and analyt- 



ical investigations of many previous researchers (Dagan 



Meier et al. 


1998: 


Sanchez- Vila et al. 


1999 


Wu et al. 


Straface et al. 


2007 


Trinchero et al. 


2008) 


, a review 



, |1982| \Ramey 




1982 




Butler 




1990 


2005| \Knudby and Carrera 




2006 



presented in the background section. The proposed approach for accounting aquifer hetero- 
geneity on Theis parameters will be presented in the methodology section. The approach 
will be demonstrated on the LANL dataset in the results section. 



2 Background 



It has been recognized that aquifer property inferences based on the Theis method and 



Cooper- Jacob method differ (Raraey, 1982). This is due to the fact that the inference meth- 



ods emphasize properties in different regions of the aquifer. The Theis method considers the 
entire drawdown curve, often leading to an emphasis on the interval of greatest curvature 



located during early times. As indicated by Butler (1990), drawdown at early times can be 



affected by many factors, including local heterogeneities near the pumping well and well skin 
and pumping storage affects, creating greater variability in Theis method inferences. The 
Cooper- Jacob method ignores early times, providing information on the properties of the 
aquifer within a ring formed by the outward moving front of the cone of depression during 
the time interval under consideration. At late time, when the Cooper-Jacob approximation 



is valid, the region included in this ring can be large. Butler (1990) demonstrates that the 



difference between inferences obtained from the Theis and Cooper- Jacob methods depend on 
the level of aquifer heterogeneity and the distance between the pumping well and the obser- 
vation location. The inferences become more similar as the level of heterogeneity decreases 
and the distance increases. 

Meier et al.\ ( 1998[ ) explore the use of the Cooper- Jacob approximation to infer effective 



transmissivity (T e ff) from the estimated transmissivity parameter T and provide indications 
of hydraulic connectivity by evaluating the estimated storativity parameter S in heteroge- 
neous aquifers. Consistent with theoretical findings of Butler ( 1990[ ), Meier et al. (1998) 
present cases where field data demonstrate that although small-scale (point) estimates of 
transmissivity T are highly variable, values of T obtained from the Cooper- Jacob method 
are relatively constant. Furthermore, Meier et al. (1998) demonstrate that corresponding 



values of S are typically highly variable, even though the storativity in the field is believed 



to be relatively constant. Meier et al. (1998) investigate this phenomena performing numeri- 



cal experiments with heterogeneous transmissivity fields and homogeneous storativity fields, 
producing similar values for T and variable values for S consistent with field cases. 

The reason for this paradoxical result can be explained by examining the equation for 
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estimating T from the Cooper- Jacob method; T = 2.3Q/AtcI, where Q is a constant pumping 
rate and I is the slope of the late-time drawdown with respect to the (base 10) log of time 
(i.e. I = (s2 — s\)/(log t 2 — log t\), where s is drawdown and t is time). This equation 
demonstrates that T is dependent on the rate of drawdown decline, which is dependent on 
the choice of t 2 and t\. Considering only the late-time drawdown where the data approximate 
a straight line with respect to log time, in accordance with the Cooper- Jacob method, means 
that the T will approximate T e jf described by the rate of drawdown after the drawdown cone 
of depression has passed the monitoring well. Storativity estimates using the Cooper- Jacob 
method (defined as S = 2.25Tto/r 2 , where r is the distance from the pumping well to the 
observation point and to represents the time-axis intercept of a line drawn through the late- 
time drawdown), on the other hand, are dependent on the variability of T between the 
pumping well and the front of the cone of depression. Although the heterogeneity between 
the pumping and monitoring well does not affect the slope of the late-time drawdown used 
to determine T, it can affect S as the time-axis intercept (to) is dependent on the arrival of 
the cone of depression at the monitoring well. If a region of high T connects the monitoring 
well and the pumping well, the value of to will be relatively small and vice-versa. As noted 
by Sanchez- Vila et al. (1999), the Cooper- Jacob method interprets an early/late arrival of a 
drawdown cone of depression as low/high storativity. This explains the high variability of S 
in the presence of T heterogeneity between the pumping well and the monitoring well, even 
in cases where S is known to be co nstant. 

Research by Meier et al. (1998) demonstrate from a numerical analysis that T estimated 
from a simulated pumping test (radial flow) is close to T e ff for parallel flow for an area of 
influence for multilognormal stationary (geostatistically homogeneous) T fields (the S field 
is assumed uniform in all cases). While Meier et al. (1998) also demonstrated that this 



can b e true for nonmultigaussian fiel ds, this i s not necessarily true in general \Sanchez-Vila 
Similar to findings by 



et al. 



1996). 



Butler 



(1988), who demonstrated that S depends on 



transmissivi ties between the pumping well and the front of the cone of depression, Meier 



et al. 



(1998) find that S depends on transmissivities between and nearby the well and the 



observation point. 



Sanchez-Vila et al. (1999) verify these conclusions using an analytical approximation to 
the groundwater flow equation. They demonstrate analytically that T is independent of 
spatial location. They also demonstrate that storativity estimates will provide an indication 
of the local deviations of T from its large-scale geometric mean (denoted as Tq) representing 
the equivalent geostatistically homogeneous T field. If T in a specific location is less than Tq, 
S will be larger than the true value of S and vice-versa. They also show that the geometric 
mean of S values is an unbiased estimator of S. 



Knudby and Carrera (2006) demonstrate that Cooper- Jacob estim ates of diffusivity (D = 

J2008b 



Trinchero et al. 



T/S) correlate well with indicators of flow and transport connectivity, 
demonstrate that estimated effective porosity (a transport connectivity indicator) depends 
on a weighted function of actual transmissivity and the interpreted Cooper- Jacob storativity 
along the flow line. 

In contrast to Meier et al. (1998), Sanchez- Vila et al. (1999), Knudby and Carrera (2006), 
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and Trinchero et al. (2008 ), Wu et al. (2005 ) explore the effect of the homogeneity assumption 



of the Theis solution on parameter estimates for the entire drawdown curve (including early 
and late time data) for cases with heterogeneous T and S fields. Conceptualizing T and S 



as spatial stochastic processes in the equation of flow, Wu et al. (2005) derive the mean flow 
equation of a heterogeneous confined aquifer as 



T eff V 2 (h) = S, 



d(h) 

dt 



X 



where angle brackets indicate ensemble mean, t is time, T e ff is the effective transmissivity 
defined as 



T, 



eff 



T + 



(T'Vti) 



(2) 



and S e ff is the effective storativity defined as 



S, 



eff 



s + 



C'dh' 
D dt 

9(h) 
dt 



(3) 



where the over bar and prime denote the spatial mean and perturbation of the variable, 
respectively. T e ff and S e ff are denoted as effective parameters as they will produce the 
ensemble mean head (h) as a convergent average for a set of realizations of heterogeneity 
based on the stochastic parameters T = T + T' and S = S + S' . As indicated by Wu et al. 
(2005), in order for the ensemble mean head (h) to equal the spatially averaged head h 



of a single realization of heterogeneity, the field must contain an adequate sampling of the 
heterogeneity (i.e. the field must be ergodic). As traditional pumping tests typically estimate 
T e ff and S e ff based on one or a small number of point estimates of head, which will not 
equal the spatially averaged head in general, T and S will not provide estimates of effective 
properties in an ensemble sense in general. 



Wu et al. (2005) performed numerical experiments using synthetic aquifers with multi- 
lognormal heterogeneous T and S fields. They observe that at early time, T estimates at 



different locations are highly variable, while, similar to the findings of Meier et al. (1998) 



and Sanchez-Vila et al. (1999), at large times (when the Cooper- Jacob approximation is 



valid) values of T converge to a value close to Tq as the cone of depression expands for the 
multilognormal fields considered. As the considered transmissivity field is multilognormal, 
Tq = T e ff. In cases where the transmissivity is nonmultigaussian, the significance of T is 



less certain { Sanchez-Vila et al. L 1999 ) , howe ver, we assume that it is a good first estimate of 
T e ff. In the same analysis, Wu et al. (2005) demonstrated that values of 5* do not converge 
to a single value, but stabilize relatively quickly to values predominantly dependent on the 
heterogeneity between the pumping well and the given monitoring location. 



Similarities to the numerical results of Wu et al. (2005) can be seen in the analytical 



investigation by Dagan (1982), who utilized a perturbation expansion approach to explore 
the temporal behavior of K e ff = — (q)/V(/i), (T e /f in 2D), where q is discharge. He derived 
an approximate relation describing the temporal behavior of T e ff for the idealized case of 
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sufficiently small transmissivity variance and average head gradient slowly varying spatially 
and temporally in a stationary random field as 



Teff(t) 



1 + -a 
2 



2 b 2 



it) 



(4) 



where \iy is the mean and a\ is the variance of Y — ln(T) and b 2 (t) is a function describing 
the temporal dependency of T e ff based on the aquifer heterogeneity in 2D, equal to unity 
for t = and tending to zero as t — > oo. Recognizing that the limiting cases for equation Q 
are first-order approximations of the arithmetic mean transmissivity Ta (t = 0) and T e ff 
{t — > oo), b 2 if) can be expressed as 



b 2 (t) 



Teff(t)-T ( 



eff,c 



Ta - T eff)C 



(5) 



where T e ^ c is the late-time convergent T e fj. Equation (JsJ) indicates that b 2 (t) describes the 
temporal decline of T e ff from Ta to T e // jC . 

In contrast to the four field cases discussed by Meier et al. (1998) (i.e. Grimsel test 
site, Switzerland {Frick 1992); El Cabril site, Spain (Bureau de Recherches Geologiques 



et Minieres); Horkheimer Insel site, Germany (Schad and Teutsch 



Air Force Base, U.S.A. (Herweijer and Young, 1991)), Straface et al. 



1994); and Columbus 



(2007) observe a lack 



of similar slope for drawdown vs log time at late times from pumping tests near Montalto 
Uffugo Scalo, Italy, indicating that the Cooper- Jacob straight-line approximation for late- 
time drawdown will not be valid in all cases. Based on their analysis of these pumping tests, 



Straface et al. (2007) question the validity of using traditional pumping tests to estimate 



meaningful hydrogeological parameters, but do suggest that these results can provide quick 
inexpensive first estimates. Furthermore, they suggest that these first estimates can be useful 
as starting parameters for a tomographic inversion of the same dataset. 



Harp and Vesselinov (2010) demonstrate an approach to identify and decompose the 
pressure influences at a monitoring location using the Theis solution. Their approach is 
demonstrated on the same dataset as in the current research. As the objective of the research 
in Harp and Vesselinov (2010) is the decomposition of pressure influences, attempts are not 
made to account for early time behavior of the Theis solution parameters, and constant 
and distinct values are applied to pumping/monitoring well pairs. Therefore, the parameter 
estimates are not considered representative of the aquifer properties of the aquifer, but are 
interpreted parameters characterizing the hydraulic response at a monitoring location due to 
pumping a single well. These interpreted parameters are analogous to parameter estimates 
that would be obtained from a conventional pumping test analysis. 

The current research presents an approach to account for Theis parameter behavior to 
infer aquifer properties considering the extensive body of research presented above. While 
the current approach is demonstrated on the long-term dataset from the LANL site, pro- 



viding the decomposition of pressure influences similar to the approach presented in Harp 



and Vesselinov (2010), the current approach could also be applied to a conventional pump- 



ing test to more appropriately account for the behavior of the Theis solution parameters. 
Furthermore, this could be particularly useful to obtain late time hydrogeologic inferences 
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from conventional pumping tests that were not conducted for a sufficient length of time to 
establish quasi-steady state conditions. 



3 Methodology 

The Theis solution of the flow equation in homogeneous media (TV 2 /i = Sdh/dt) is defined 
as 

8 (f) = -®-W{u) = -Q-W ( r ^f) , (6) 
py ' 4vrT v ; 4ttT \4Tt /' w 

where s p (t) is the predicted pumping drawdown at time t since the pumping commenced (i.e. 
h(t) —h(0)), Q is the pumping rate, T is the transmissivity, W(u) is the negative exponential 
integral (/^° e~ y /ydy) referred to as the well function, u is a dimensionless variable, r is radial 
distance from the pumping well, and S is the storativity. Multiple pumping wells and variable 
rate pumping periods can be included in the Theis solution by employing the principle of 



superposition {Freeze and Cherry, 1979, page 327) as 



s 



{t)=j:Y: Qi \ ^ wi ^s X (7) 



=1 j: 



^ AnT \4T(t — t 



where N is the number of pumping wells (sources), M, is the number of pumping periods 
(i.e. the number of pumping rate changes) for pumping well i, Qij is the pumping rate 
of the ith well during the jth pumping period, and tg i:j is the time when the pumping 
rate changed at the ith well to the jth pumping period. The drawdown calculated by 
equation ([7| represents the cumulative influence at a monitoring location of the N pumping 
wells at distances r iy i = 1, . . . , N from the monitoring location. 

Equations ^ and ^ are only valid under the assumption of homogeneity. If a system 
is homogeneous, then T and S in equations ^ and ^ will be equivalent to T e ff and S e ff, 
respectively. If the system is heterogeneous, this will only be true in an ensemble mean sense. 
In this case, the Theis solution can be expressed as 



which is the solution to equation Q ( Wu et al, 2005), where (s p )(t) is the ensemble mean 



drawdown due to pumping and heads have been converted to drawdown. Invoking superpo- 
sition with equation Q, an ensemble mean drawdown equation analogous to equation ([7]) 
can be expressed as 

As water elevations recorded at monitoring wells in an aquifer system are merely point 
samples from a single realization of heterogeneity, and not ensemble mean values of multiple 
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realizations or spatial averages of an ergodic field, application of equations ^ and ^ are 
invalid for cross-hole interference tests. Recognizing this theoretical limitation of applying 
the Theis solution (or the Cooper- Jacob approximation) to data from heterogeneous aquifers 
to infer effective properties, researchers have investigated what information is contained in 



the hydrogeologic parameter estimates ( Meier et al. , 1998 ; Sanchez- Vila et al. , 1999 Wu 



et al, 2005 Knudby and Carrera, 2006 Trinchero et al. 2008). We propose that although 



the Theis solution parameters will not provide precise representations of hydrogeological 
properties, the analytical framework of the Theis solution can provide initial estimates of 
the effective transmissivity and indications of connectivity. 

Recognizing that a dataset containing drawdown outside of the Cooper- Jacob domain 
will require consideration of the behavior of parameter estimates at early times (as the front 
of the cone of depression is at short radial distance), we approximate the Theis solution, 
defining the estimated pumping drawdown s p (t) as 



N 



EE g 

i=i j=i 



Qij 



-w 



m * {t 




(10) 



where and Si are time dependent functions describing the variation in interpreted trans- 
missivities and storativities as the cone of depression propagates outward from the pumping 
well. In order to provide a general functional form with the intent to capture the temporal 
dependence of T and S for a broad range of heterogeneities and pumping well factors, T, 
and Si are defined using an exponential functional form as 



Ti(t) 



■eff 



CT > 0, 



'ID 



and 



Si{t) 



q ,cs,i/{t~tQ l ,) 



(12) 



where T e // provides the late-time convergent estimate for T e ff, S a> i provides late-time con- 
vergent indications of connectivity between the ith pumping well and the monitoring location 
(Knudby and Carrera, 2006), and c^i and c$i are constants describing the time dependent 



slope of the transmissivity and storativity parameters, respectively, associated with the ith 
pumping well. Since in most cases, the statistical nature of the heterogeneity will not be 
known with certainty, this ad hoc functional form is assumed reasonable. In idealized sce- 
narios with known correlation structure, it may be possible to derive these relationships 
in an ensemble mean sense. For example, Dagan (1982) derives an analytical relationship 
describing the temporal behavi or of T e f f f or an exponential autocorrelation. 
Based on 



Dagan (1982) and 



Wu et al. 



(2005), we constrain ct > 0, indicating that T{t) 



values from early time portions of drawdown data are expected to be higher than late time 
convergent values. This may be explained by the early-time negative correlation between 



head and transmissivity ( Wu et al. , 2005) and/or the existence of high conductivity inter- well 
pathways as described by Herweijer (1996). Other possible explanations for time-dependent 
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hydrogeologic parameters are well-bore storage and leakage effects known to exist at the site 
flMc£ml|2005l|2006a|b) . 



Substituting equations (Jllp and (JT2j) into equation ( J10 ) produces 

Qi,j Qij—1 t t r ( '"i ^a,i^ 



N Mi 



i=l j=l 47TJ e //e ^ 



(13) 



In order to account for temporal trends identified in a previous study (Harp and Vesseli- 



nov 



2010), we include an additional drawdown term s t (t) as 



s t {t) = (t-U)xm (14) 

where t% is the time at the beginning of the pumping record and m is a constant defining the 
linear increase in drawdown not attributable to pumping. 

As the calibration targets in the model inversions presented here are water elevations as 
opposed to drawdowns, we define the predicted water elevation h(t) at time t as 

k{t) = h - s p (t) - s t (t) (15) 

where h Q = h(0) and is defined as the initial predicted water elevation at the monitoring well. 
As defined by the Theis solution, h is the head at the time that a perturbation commences. 
As pumping of the regional aquifer began at the LANL site over 50 years ago, it is reasonable 
to assume that the influence of the earlier pumping has propagated through the system 
and/or dissipated. However, more recent pumping rate changes preceding pressure transient 
records at monitoring locations need to be considered. In order to account for residual effects 
of pumping prior to monitoring, simulations are started in advance of pressure transient 
records including earlier pumping records. Therefore, h Q is not a measured quantity, but 
predicted at the beginning of the simulations. 



Model calibration is performed using a Levenberg-Marquardt approach (Levenberg, 1944 



Marquardt, 1963) where the objective function can be defined as 



h t (t 3 )f 



(16) 



i=l j=l 

where m is the number of monitoring wells considered, is the number of head observations 
for the zth monitoring well, and hi(tj) are the head observations for the ith monitoring well 
included as calibration targets where j is an observation time index. 

The simulation of the drawdowns is performed using the WELLS code (available upon 
request at http://www.ees.lanl.gov/staff/monty/), which implements equation (13). 



calibration is performed using PEST ( \Doherty 2004) 



The 



4 Site Data 

The regional aquifer beneath the LANL site is a complex stratified hydrogeologic structure 
which includes unconfined zones (under phreatic conditions near the regional water table) and 
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confined zones (deeper zones) ( Vesselinov 2004af b). The three monitoring wells considered 
in this analysis are screened near the top of the aquifer with an average screen length of 11 
meters. The water-supply wells partially penetrate the regional aquifer with screens that 
also begin near the top of the aquifer, but penetrate deeper with an average screen length 
of 464 meters. Nevertheless, field tests demonstrate that most of the groundwater supply 
is produced from a relatively narrow section of the regional aquifer that is about 200-300 
m below the regional water table {Los Alamos National Laboratory, 2008a). Implicit in 



the use of the Theis solution is the assumption that groundwater flow is confined and two- 
dimensional. We assume that this is a justifiable assumption here given the small magnitude 
of observed drawdowns (less than 1 m at the monitoring wells and less than 20 m at the 
water-supply wells) and the relative long distances between supply and monitoring wells 
compared to the effective aquifer thickness (about 200-300 m). It is believed that leakage 
and vertical flow may be significant factors within the aquifer system, however, lack of 
monitoring locations within aquifers and confining units and lack of sufficient information 
to identify the locations of aquifers and confining units hinders approaches that intend to 
account for these affects (e.g. \Neuman and Witherspoon (1972)). 

Water-level fluctuations (pressure transients) are automatically monitored using pressure 
transducers. The pressure and water-supply pumping records considered here are collected 
from 3 monitoring wells (R-H, R-15 and R-28) and the 7 water-supply wells (PM-1, PM-2, 
PM-3, PM-4, PM-5, O-l, and 0-4) located within the LANL site. Figure [l] displays a map 
of the relative location of the wells. Figure [2] presents the pressure and production records 
for the monitoring wells and water-supply wells, respectively. 

The water-level observation data considered here span approximately five years, com- 
mencing on or shortly after the date of installation of pressure transducers (May 4, 2005 for 
R-ll; December 23, 2004 for R-15; February 14, 2005 for R-28), all terminating on Novem- 
ber 20, 2009. The barometric pressure fluctuations are removed using constant coefficient 



methods using 100% barometric efficiency for all monitoring wells {Los Alamos National Lab- 



oratory 2008b). Although the pressure transducers collect observations every 15 minutes, 



this dataset is reduced to single daily observations by using the earliest recorded measure- 
ment for each day. Some daily observations have been excluded due to equipment failure. 
Pumping records for all pumping wells begin on October 9, 2004 and terminate on December 
31, 2009. The pumping record precedes the water-level calibration data to account for water- 
level transients due to pumping variations before the water-level data collection commenced. 
For additional information on the site and dataset, refer to Harp and Vesselinov (2010). 

In the applied computational framework, forward model run times for predicting water 
elevations at R-ll, R-15, and R-28 for approximately four years (from October 8, 2004 
to November 18, 2008) are each approximately 3 seconds on a 3.0 GHz Intel processor. 
Inversions initiated with uniform initial parameter values require approximately 600 model 
runs and, using a single processor, are performed for approximately 1 hour and 40 minutes. 



10 



0-1 

★ 



0-4 
★ 



PM-5 
★ 



PM-3 
★ 



PM-1 
★ 



PM-4 
★ 



PM-2 
★ 



Okm 0.5 km 1km 1.5 km 2 km 



Figure 1: Map of monitoring wells (circles) and water-supply wells (stars) included in the 
analysis. Locations of newly completed and planned monitoring wells are indicated by open 
diamonds. 
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Figure 2: Water elevations at monitoring wells and production records for water-supply 
wells. 
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Figure 3: Calibrated heads for the exponential (red) and constant (black) cases. The ob- 
served heads are presented in gray. 



5 Results and Discussion 

Results of calibrations using temporally varying parameters and constant parameters are 
presented below. In the exponential case, a single value of T e ff is applied to all pump- 
ing/monitoring well pairs. Similarly, in the constant case, a single value of T is applied to all 
pumping/ monitoring well pairs. Distinct values are allowed for S a and S in the exponential 
and constant cases, respectively. Pumping influences (wells) that the calibration is unable 
to fingerprint at the monitoring well result in unrealistic parameter values that effectively 
eliminate the influence of the pumping well (i.e. high T and S). As these parameter val- 
ues are not physically meaningful beyond identifying a lack of influence from the associated 
pumping well, they are not presented below. Therefore, the omission of a pumping well 
below indicates a lack of identifiable influence at a monitoring well by the inversion. 

Figures [3] presents the calibrated drawdowns from the water-supply wells for monitoring 
wells R-ll, R-15, and R-28 for the exponential and constant cases. It is apparent that 
the exponential case reduces the mismatch in drawdowns for all three wells, with the most 
significant improvements in R-15, the well with the worst match for the constant case. 

Figure [4] presents the estimated transmissivity functions (equation 
and R-28. The functions are plotted up to around five years to include parameter values 
used in the model runs. It is apparent that for all three monitoring wells, T(t) converges 
towards a single value, T e ff = 10 3 07 m 2 /d (1170 m 2 /d), as constrained by the inversion. The 
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Figure 4: Estimated transmissivity functions for the exponential case. The convergent value 
of T e ff and T for the constant case are indicated. 
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Figure 5: Estimated storativity functions for the exponential case. 
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interpreted transmissivity for the constant case (T = 10 3 ' 27 m 2 /d (1860 m 2 /d)) is indicated 
in the figure. It is apparent that this value is fitted to an average value of T(t) for the 
exponential case. This indicates that an overestimate of T e ff will be obtained using constant 
parameters. 

Figure [5] presents the estimated storativity functions (equation (12)). It is apparent that 
in general the storativity functions converge quickly to distinct values in accordance with 
previous research ( Wu et al, 2005 Straface et al, 2007), providing indications of inter- well 
connectivity. Physically unrealistic values of storativity are allowed as S a is recognized as 
a flow connectivity indicator, and does not represent aquifer storativity in an effective or 
equivalent sense. 

Table [T] presents the estimated parameters associated with the transmissivity and stora- 
tivity functions plotted in Figures [4] and [5] for the exponential case and the parameters for 
the constant case. As constrained in the inversion, all transmissivities converge to a single 
value for both the exponential and constant cases. For the exponential case, this value can 
be considered a first estimate of T e ff. The larger value obtained for the constant case indi- 
cates that the calibration has fitted the parameter within the early time variability, thereby 
overestimating T e ff. 

Values of S a indicate the level of connectivity between the monitoring and pumping well. 
Large/small values of S a indicate a region of relatively low/high inter-well transmissivity. It 



is apparent from Figure |5| and Table [1| that the trends in S a can be grouped by the associ- 
ated pumping well. For instance, convergent values of S a decrease (inter-well connectivity 
increases) from PM-2 to PM-3 to PM-4 and PM-5. In general, similar trends are apparent 
for S in the constant case as well. However, in the constant case, values for PM-2 are farther 
from physically realistic values of storativity. 

A decomposition of the pressure influences from the pumping wells at the monitoring wells 
also resulted from this research. These results are similar to the decomposition analysis of 
this dataset presented in Harp and Vesselinov (2010) and therefore are not presented here. 
For instance, the same pumping wells are identified to influence drawdown at the monitoring 
wells and a lack of a linear temporal trend is identified for R-15 in both cases. 



6 Conclusions 



This paper demonstrates an approach to obtain late-time aquifer property inferences con- 
sistent with the Cooper- Jacob method from transient datasets collected in heterogeneous 
aquifers. Such datasets are commonly available from municipal water-supply networks. The 
utilization of these existing datasets eliminates the expense and coordination necessary to 
perform dedicated pumping tests at a site. The methodology is motivated by analytical 
investigations by Dagan (1982), numerical experiments by Wu et al. (2005), and analysis 



of field-collected hydrographs by Straface et al. (2007). The hydrogeologic inferences are 



evaluated based on a large body of research into the meaning of late-time aquifer prop- 



erty inferences ( 
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-116.9 
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-3.19 
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-4.7 








PM-2 






287.9 
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-433.2 
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33.3 


0.19 


0.27 


-19.4 
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0.02 




PM-4 






29.4 


0.08 


0.13 


-21.3 







Table 1: Parameter estimates from calibrations using exponential functions (Exp.) and con- 
stant (Const.) parameters. Transmissivity parameters are presented in units of log(m 2 /d). 
m is the linear temporal trend parameter. 



Utilizing this approach on a dataset from the LANL site has indicated that adequate 
water-level calibrations can be achieved within the constraints of the inversion: a single 
value of T e ff is applied to all pumping/monitoring well pairs; T(t) decreases towards a 
constant value; S a is allowed to take distinct values and is allowed to increase or decrease 
towards convergent values. T e ff provides an initial estimate of the effective transmissivity at 
the support scale characterized by the distances between the pumping and observation wells 



( Neuman 


1990 
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2003 


)• I 
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(1999 
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of inter-well connectivity, indicating the degree in w 
are hydraulically connected. 



n acco rdance with Meier et al. (1998), 
( 2006 ) , S a is recognized as an indicator 
rich pumping and monitoring well pairs 
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